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1 Introduction 

(N 

\Q \ The importance of the strange quark mass, as a fundamental parameter of 

^ \ the Standard Model (SM) and as an input to many interesting quantities, 

QQ ' has been highlighted in many reviews, eg in Ref [jll. A first principles 

O ■ calculation of is possible in lattice QCD but to date there has been a 

. rather large spread in values from lattice calculations. This review aims 

to clarify the situation by explaining the particular systematic errors and 
Q^l their effects and illustrating the emerging consenus. 

Q^' In addition, a discussion of the strange quark mass is timely given the 

^ ■ recent results from KTeV and NA48 for e'/e which firmly establish 

direct CP-violation in the SM and when combined with previous measure- 
ments give a world average e'/e = (21.2 ± 2.8) x 10~^. This is in stark 
^ \ disagreement with the theoretical predictions which favour a low e'/e [Q. 

' Although in principle e'/e does not depend directly on rris in practice it 

has been an input in current phenomenological analyses. This dependence 
arises because the matrix elements of the gluonic, (Qe)^ and electroweak, 
{Qs)2^ penguin operators]^ are of the form (7r7r|Qj|i^) and final state interac- 
tions make them notoriously difficult to calculate directly. They have been, 
therefore, parameterised in terms of bag parameters, the strange quark 
mass, rris and the top quark mass, mj, as discussed in detail in Ref. 0. A 
recent review of lattice calculations of the matrix elements is in Ref. 0. In 
this talk I will focus on some recent and careful lattice determinations of 
rris, illustrating the reasons for the large spread in earlier results. 



*Talk presented at KAON'99, University of Chicago, June 1999 
^keeping only the numerically dominant contributions for simplicity 
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2 The stange quark mass from lattice QCD 



In lattice QCD, is determined in two ways, each of which rehes on 
calculations of experimentally measured quantities to fix the lattice bare 
coupling and quark masses. The IP- IS Charmonium splitting, Mp and 
ro are some of the parameters typically chosen to fix the inverse lattice 
spacing, a~^. To determine either Mk or is used. It is an artefact 
of the quenched approximation that depends on the choice of input 
parameters, so that some of the spread in answers from lattice QCD can 
be attributed to different choices here. Naturally, some quantities are bet- 
ter choices than others being less sensitive to quenching or having smaller 
systematic errors. 

The quark mass can be determined from hadron spectroscopy, using 
chiral perturbation theory to match a lattice calculation of Mk (or M^) to 
its experimental value with 

Ks = Bps ^ ^ +■■■ or My = Ay + Bv ^ ^ + . . . (1) 

This is the hadron spectrum or vector Ward identity (VWI) method. 

Alternatively, the axial Ward identity (AWI): d^A^{x) = {mi + mj)P{x) 
imposed at quark masses to correspond to either the experimentally mea- 
sured Mk or determines m^. 

The lattice bare masses and matrix elements are related to their con- 
tinuum counterparts, in say the MS scheme, by the renormalisation coef- 
ficients, Zg or p), calculated perturbatively or nonperturbatively. 



^ Zp[ijL,ma) {F{x)J [}})) 



rris has been calculated in all three lattice fermion formalisms: Wilson, 
staggered and domain wall. Although the domain wall fermion results are 
extremely interesting, since this approach has the good flavour structure 
of Wilson fermions while preserving chiral symmetry, the results for 
are still preliminary so I will focus on results with Wilson and staggered 
fermions. A description of the domain wall formalism and results can be 
found in Ref. p. 

Comparing results from these different methods provides a nice check 
of lattice calculations. 



3 Main uncertainties in the calculation 

The difference in early lattice results can be understood in terms of the 
treatment of systematic uncertainties in these particular calculations. The 
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largest of these are discretisation errors, calculation of renormalisation co- 
efficients and the quenched approximation. 

1. Discretisation Errors : The Wilson action has discretisation errors of 
0(a), so for a reliable result one needs fine lattices and a continuum extrap- 
olation, a ^ 0. See Figure |l] for the CP-PACS collaboration's quenched 
Wilson results |p. The Sheikholeslami-Wohlert (SW) clover action in- 
cludes a term ~ csw^^iJ.uF^y^ and discretisation errors start at O^aso), 
when csw is determined perturbatively. The remaining a-dependence must 
be removed by continuum extrapolation, but the slope of the extrapola- 
tion is milder A nonperturbative determination of csw [0 gives an 
(9(a)-improved action, which should futher reduce the lattice spacing de- 
pendence. Recent results from the APE, ALPHA/UKQCD and QCDSF 
collaborations use this approach [|1^, |15|, [l^. The latter two groups in- 
clude continuum extrapolations and find significant a-dependence (~ 15% 
between the finest lattice and a = as found by ALPHA/UKQCD). In the 
case of the more commonly used VWI approach the slope of the extrapo- 
lation in a is positive and therefore, at finite lattice spacing is too high, 
even with improvement. 

The staggered fermion action is (9(a)-improved so the lattice spacing 
dependence should be mild. 

2. Renormalisation coefficients : and -Z^(a,p) can be determined per- 
turbatively or nonperturbatively. A nonperturbative calculation is prefer- 
able as it removes any perturbative ambiguity. This was pioneered by the 



APE and ALPHA groups pO, 11 



For Wilson fermions perturbative corrections are smaller and therefore 
more reliable in the VWI approach (ie. for Zs) than in the AWI approach. 
In Ref. |jl2| the difference between nonperturbative results and boosted 
perturbation theory is ~ 10% for Zs and ~ 30% for Zp at ~ 2.6 
GeV. For staggered fermions the perturbative coefficients are large and 
positive so the results are unreliable and nonperturbative renormalisation 
is essential. The perturbative staggered results are therefore too low and 
this effect combined with the too high values of m<j from Wilson results at 
finite lattice spacing explain much of the spread in lattice results. 

3. Quenching : Most calculations are done in the quenched approxi- 
mation - neglecting internal quark loops - as a computational expedient. 
An estimate of this approximation, based on phenomenological arguments, 
was made in 0. The authors estimated that unquenching lowers by 
pa 20 — 40%. They also argued that Mk rather than is a better choice 
of input parameter since it is less sensitive to quenching. Unquenched cal- 
culations by CP-PACS have shown that these estimates were of the correct 
size and sign [ P^ . 

A number of clear trends are therefore identified: 
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• There is significant a-dependence in the Wilson action resuhs which 
raises iris at finite lattice spacing. Although this is milder for the 
improved actions it is still present, as pointed out in Refs. |T^. 

• Using perturbative improvement, the VWI and AWI methods differ 
at finite lattice spacing but agree after continuum extrapolation. This 
indicates the methods have discretisation errors larger than the per- 
turbative uncertainty. Nonperturbative renormalistion has a larger 
effect on AWI results, bringing them into agreement with VWI re- 
sults at finite lattice spacing. However, discretisation errors remain 
a significant uncertainty and without a continuum extrapolation lead 
to an overestimate of m^. 

• Perturbative renormlisation of staggered fermions results in an un- 
derestimate of nis- Nonperturbative renormalisation is essential. 

• A lower value of rris is expected from an unquenched calculation. 



4 Recent results for iris 

The systematic uncertainties in the lattice determination of are now 
well understood. Some recent results which I believe provide a definitive 
value of Tils in quenched QCD and an unquenched result are now discussed. 



4.1 Quenched results 



Table compares a number of recent calculations of m^. The JLQCD, AL- 
PHA/UKQCD and QCDSF groups have removed all uncertainties within 
the quenched approximation. JLQCD use staggered fermions and nonper- 
turbative renormalisation They observe mild a-dependence, as ex- 
pected and take the continuum limit. The effect of nonperturbative renor- 
malisation is considerable, again as expected: ~ +18% when compared to 
the perturbative result. 

The ALPHA/UKQCD collaborations Jll and QCDSF use a non- 



perturbatively improved SW action and renormalisation and include a con- 
tinuum extrapolation. This explains the difference between their results 
and that of the APE group (which has not been extrapolated to a = 0). In- 
terestingly, ALPHA/UKQCD, QCDSF and the Fermilab g and LANL 
results for are in agreement. The difference in analyses is nonperturba- 
tive versus perturbative renormalisation, indicating that the perturbative 
result for the VWI method is reliable (for Wilson fermions). 
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GROUP 


# lattice spacmgs a — > U (/OeV) 


APE '98 [g 


2 NO 121(13) 


FNAL '96 
LANL '96 


3 YES 95(16) 
3 YES 100(21)(10) 


JLQCD '99 PI 
QCDSF '99 llj 
ALPHA/UKQCD '99 |15| 


4 YES 106(7.1) 

3 YES 105(4) 

4 YES 97(4) 



Table 1: Quenched lattice results, the APE result is obtained at 
a-^=2.6GeV 

4.2 An unquenched result 

There are a number of new preliminary unquenched calculations of rUg [0 
however, CP-PACS have recently completed their analysis [l^, shown in 
Figure |1|, so I will concentrate on this. Since unquenching requires a huge 
increase in computing time it is prudent to use coarser (less time consum- 
ing) lattice spacings. This in turn requires improved actions to control 
the discretisation effects. CP-PACS use a perturbatively improved quark 
and gluon action and extrapolate to the continuum limit. The perturbative 
renormalisation is reliable with a remaining perturbative error of (9(a^), for 
Wilson fermions. The final result is m*-^'^(2GeV) = 84(7)MeV. Although 



strange quark mass, K input 



strange quark mass, input 




^ 140 

> 
o 

a 120 

cj 
3 

e" 100 















>--'"' 













0.4 0.6 0.8 

a [GeV"'] 

♦ RC, VWI with sea quark 

• RC, VWI with PQ 
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a [GeV"'] 

♦ RC, VWI with sea quark 

• RC, VWI with PQ 
■ RC, AWI 

o quenched PW, VWI 
□ quenched PW, AWI 



Figure 1: RC is the RG-improved gluon and Clover action (with two definitions 
of quark mass in the chiral limit), PW are the quenched Wilson results already 
discussed. 

this result disagrees with bounds derived from the positivity of the spectral 



function ||T8| it remains unclear at what scale, fx, perturbative QCD and 
thus the bound itself becomes reliable. CP-PACS conclude that unquench- 
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ing lowers rus - compare the filled and open symbols in Figure ^ As in the 
quenched case the VWI and AWI methods differ at finite lattice spacing 
but extrapolate to the same result - compare the ^ symbols. Finally, 
the strange quark mass obtained from the K and mesons yields consistent 
continuum values in full QCD: 84(7) MeV and 87(11) MeV respectively. 



5 Conclusions 

There has been much progress this year in lattice calculations of m^. Cur- 
rent computing power and theoretical understanding are sufficient to deter- 
mine rris to great precision. A calculation removing all uncertainties would 
include unquenched simulations, a continuum extrapolation and nonper- 
turbative renormalisation and can be done in the short term. Simulations 
at n/ = 2 and 4 with an interpolation to nj = 3 are also a possibility. 

Finally, I look at the implications for e'/e from current theoretical cal- 
culations given the recent lattice calculations of m^. The dependence is 
shown in Figure from the analytic expression 



e'/e = ImXt 



Co + (ceS?/'^ + csBs 



(3/2) \ 



K 



(2) 



and input from lattice calculations for Bi [T^. The values of other SM 
parameters are from Ref. |2^. The lines represent the effect of varying 
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Figure 2: e'/e as a function of from Equation. 



the bag parameters and/or the Wilson coefficients and the band is the un- 
quenched iTis from CP-PACS, run to rric. Further reducing the uncertainty 
on iTLs is more straightforward than for the Bi and can constrain theoretical 
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calculations of e'/e. Clearly the lower values of give higher e'/e values, 
in better agreement with experiment! 
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